'3>T\C 


The  optical  properties  of  aerosols 


Final  Technical  Report 
by 

F.  Borghese(  Q  o  ^  VYVeiSvrcx 
August  1998 


.  1  b6-A-^ 


United  States  Army 

EUROPEAN  RESEARCH  OFFICE  OF  THE  U.  S.  ARMY 
London,  England 


Contract  number  N68171-97-M-5546 
Contractor:  Centro  Siciliano  per  le  Ricerche  Atmosferiche 
'  e  di  Fisisca  dell’Ambiente 


Approved  for  public  release;  distribution  unlimited. 


r 


The  optical  properties  of  aerosols 
Final  Technical  Report  on  Contract  N68171-97-M-5546 
August  1998 

The  research  performed  under  the  present  contract  has  been  addressed  to  the 
following  topics; 

1.  Optical  properties  of  aggregated  spheres  deposited  on  a  dielectric  surface 
with  a  random  distribution  of  their  orientations. 

2.  Determination  of  atmospheric  pollution  through  spectral  analysis  of  starlight. 
The  reserach  mentioned  above  has  been  reported  in  a  number  of  papers  that 
were  already  printed  or  are  at  present  in  the  process  of  printing. 

At  the  same  time,  our  reserach  called  our  attentiuon  on  a  number  of  mathemati¬ 
cal  problems  thatn  if  not  properly  faqed,  may  frustrate  any  attempt  to  get  relia¬ 
ble  numerical  resuts. 

The  main  findings  that' are  reported  in  our  papers  will  now  be  briefly'summar- 
ized,  and  the  lines  of  future  progress  are  sketched.  -  -  - 

1.  Optical  properties  of  aggregated  spheres  deposited  on  a  dielectric  surface 
with  a  random  distribution  of  their  orientations. 

This  research  is  meant  to  extend  to  the  case  of  aggregated  spheres  the  theory 
that  we  already  developed  to  deal  with  the  extinction  spectra  from  spheres  de¬ 
posited  on  a  dielectric  substrate.  The  need  for  such  an  extension  is  easily  under¬ 
stood  on  account  that  spherical  particles  tend  to  aggregate  to  form  compound 
scatterers  whose  otical  properties  are  known  to  be  rather  different  from  the  the 
properties  of  the  component  spheres.  Of  course,  the  motivation  of  our  study  is 
to  supply  a  reliable  method  to  check  the  cleannes  of  surfaces. 

Our  starting  point  has  been  the  theory  that  we  already  developed  to  deal  with 
the  spectra  from  aggregates  of  known  orientation.  Using  the  transformation 
properties  of  the  spherical  multipole  fields  under  rotation  we  were  able  to  de¬ 
fine  the  transition  matrix  for  an  aggregate  in  the  presence  of  a  dielectric  sub¬ 
strate  and  to  perform  analitically  the  required  averages  over  the  orientational 
distribution.  As  a  result  we  were  able  to  produce  the  full  scattering  patterns 
from  a  dispersion  of  binary  aggregates  deposited  on  a  dielectri  substrate  with  a 
random  distribution  of  their  orientations,  including  both  co-polarized  as  well  as 
cross-polarized  patterns.  The  interpretation  of  our  results  has  been  reported  in 
a  paper  that,  though  not  yet  printed,  has  been  accepted  for  publication  in  JOSA 
A. 

An  important  point  regards  the  comparison  of  the  results  of  our  theory  with  the 
experimental  data.  To  this  end,  the  experimental  group  at  the  Dipartimento  di 
Fisica  della  Materia  e  Tecnologie  Fisiche  Avanzate,  was  able  to  devise  an  exper¬ 
imental  setup  to  record  the  spectra  from  a  surface  sparsely  seeded  by  particles. 
Wen  the  resulting  recoreded  spectra  were  compared  with  the  results  of  our  cal¬ 
culations,  a  complete  agreement  was  found  within  the  experimental  error. 


19981006  035 


2.  Determination  of  atmospheric  pollution  through  spectral  analysis  of 
starlight. 

It  is  well  known  that 'the  optical  methods  are  a  useful  tool  for  the  determinan- 
tion  of  atmospheric  pollution.  As  a  rule  the  light  that  serves  as  a  probe  is  gener- 
•  ated  on  the  surface  of  Earth  as  in  LIDAR  techniques,  or  comes  from  satellites  in 
appropriate  orbits.  In  recent  times  our  attention  has  been  drown  on  a  third  pos¬ 
sibility,  i.  e.  on  the  spectral  analisis  of  the  light  that  comes  from  selected  stars. 

In  principle,  the  proposed  tecnique  is  rather  simple.  The  light  from  a  star  is  tak¬ 
en  form  two  observers  located  at  a  different  altitude  is  analized  and  the  result¬ 
ing  spectra  are  compared.  It  is  intuitive  that  any  extra  absorpion  line  that  might 
appear  in  the  spectra  taken  at  low  altitude  should  originate  from  atmospheric 
pollution.  In  practice  the  realization  of  this  kind  of  measurements  requires  the 
■presence  in  the  same  area,  of  an  astronomical  observatory  and  of  a  low  altitude 
laboratory.  A  situation  of  this  kind  occurs  in  the  area  of  Catania,  Italy,  where 
one  can  use  the  facilities  of  the  Astronomical  Observatory  on  Mount  Etna,  some.  • 
2900  metres  of  altitude,  as  well  as  the  Laboratory  of  the  University  of  Catania, 
that  is  located  almost  at  sea  level.  Note  that  the  horizontal  distance  between  the 
Observatory  and  the  Laboratory  does  not  exceed  10  kilometres.  Similar  situa¬ 
tions  occur  at  other  places  in  the  World,  e.  g.  in  Chile,  South  America. 

From  a  physical  point  of  view,  the  actual  determination  of  the  polluting  sub¬ 
stances  is  complicated  from  the  fact  that,  specially  in  the  infrared,  one  should 
take  into  account  the  effects  that  are  due  to  the  presence  of  particles.  The  effect 
of  the  latter  objects  can  be  discarded  from  the  observed  spectra  provided  that 
one  is  a'ble  to  calculate  the  extinction  spectra  from  a  polydispersion  of  particles. 

It  is  at  this  stage  that  the  tecniques  that  we  developed  under  the  present  and  the 
previous  Contracts  enter  into  play.  Indeed,  even  the  most  anisotropic  particles 
can  be  dealt  with  by  our  codes,  provided  that  the  spectroscopic  data  on  their 
chemical  composition  are  known.  At  present,  this  researech  is  performed  by  dr 
Giovanni  Catanzaro  as  a  part  of  his  work  to  get  the  degree  of  Dottore  di  Ricerca 
in  Fisica.  Of  course,  the  work  is  far  from  concluded  and  we  plan  to  pursue  the 
subject  in.  the  framework  of  the  next  Cotract  that  has  already  been  awarded  by 
the  ERO. 

I  cannot  conclude  this  report  withuot  mentioning  that  thanks  to  the  work  per¬ 
formed  under  the  present  Contract,  dr  Maria  Antonia  Iati,  who  collaborated  to 
the  work  of  the  last  two  yiears,  got  the  degree  of  Dottore  di  Ricerca  in  Fisica 
and  is  presently  performing  post  doctoral  research  at  the  University  of  Water¬ 
loo,  Ontario,  Canada,  under  the  direction  of  prof.  Duley. 
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Abstract 

The  preliminary  results  of  the  calculation  of  the  full  pattern  of  the  scattered  inten¬ 
sity  from  a  dispersion  of  randomly  oriented  binary  aggregates  of  identical  spheres 
deposited  on  a  dielectric  substrate  is  presented. 

In  a  recent  paper  we  presented  the  full  pattern  of  the  scattered  intensity  from  a 
single  sphere  deposited  on  the  plane  surface  that  separates  two  homogeneous  media 
of  different  refractive  index.1  The  theory  behind  our  results  is  based  on  the  expansion 
of  the  electromagnetic  field  into  a  series  of  spherical  multipole  fields  and  on  the 
imposition  of  the  boundary  conditions  across  the  plane  surface  and  across  the  surface 
of  the  sphere.  The  difficulty  of  imposing  the  boundary  conditions  across  the  plane 
surface  to  a  field  given  in  terms  of  spherical  multipole  fields  was  overcome  through 
the  application  of  a  general  reflection  rule  for  spherical  vector  multipole  fields  on  a 
plane  surface.2  According  to  this  rule,  the  reflection  of  a  spherical  multipole  field 
with  origin  at  P  originates  a  linear  combination  of  spherical  multipole  fields  that 
satisfy  the  radiation  condition  at  infinity,  with  origin  at  P\  the  mirror  image  of  P 
with  respect  to  the  plane  surface.  Using  this  rule,  that  clarifies  the  role  of  image 
theory  when  the  surface  is  notJ  perfectly  reflecting;  we  get  a  scattered  field  that 
has  the  correct  behavior  at  infinity.  Perhaps,  the  most  distinctive  feature  of  our 
approach  is  the  fact  that,  unlike  similar  approaches  to  the  same  problem,3,4  we  get  a 
non-vanishing  field  that  propagates  along  the  surface.  As  a  result  the  patterns  that 
we  presented  were  in  excellent  agreement  with  the  experimental  findings  and  with 
the  ab  initio  simulations.5,6 

The  purpose  of  this  paper  is  to  present  some  results  for  the  full  scattering  pat¬ 
tern  from  binary  aggregates  of  identical  spheres  deposited  on  a  plane  surface.  The 
extension  of  our  theory  to  the  case  of  aggregated  spheres  is  far  from  trivial  because 
one  has  to  consider  the  effect  of  the  mutual  interaction  of  the  aggregated  spheres:  in 
other  words,  we  had  to  modify  our  theory  to  account  for  the  dependent-scattering 
effects  in  the  presence  of  the  plane  surface. 

Of  course,  aggregated  spheres  form  anisotropic  objects  whose  patterns  are  ex¬ 
pected  to  have  a  noticeable  dependence  on  their  orientation  with  respect  to  the 
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incident  field.  The  field  scattered  by  a  dispersion  of  such  objects  even  when  de¬ 
posited  on  a  surface  with  a  given  distribution  of  their  orientations  may  be  rather 
different  from  the  field  scattered  by  a  dispersion  of  the  same  objects  all  oriented  alike. 
Nevertheless,  by  making  full  use  of  the  transformation  properties  of  the  spherical 
multipole  fields  under  rotation7  we  are  able  to  calculate,  through  an  analytical  av¬ 
erage,  the  field  scattered,  for  instance,  by  an  orientationally  random  dispersion  of 
aggregates  on  the  surface.8  This  has  been  achieved  by  giving  a  suitable  definition  of 
the  transition  matrix1,9  of  the  aggregate  in  the  presence  of  the  surface.  The  elements 
of  the  transition  matrix  include,  of  course,  all  the  interactions  among  the  spheres 
that  compose  the  aggregate.  The  resulting  patterns,  that  will  be  discussed  below, 
should  give  a  reliable  description  of  the  physical  situation  that  is  likely  to  be  meet, 
in  actual  experiments. 


In  the  figure  above  we  sketch  the  geometry  that  we  use  to  report  our  results 
for  the  scattering  patterns.  The  surface  of  the  substrate  is  the  x-z  plane  and  the 
direction  of  incidence  is  defined  by  the  polar  angles  i?  j  and  (pj .  In  particular  for  the 
results  reported  in  this  paper  we  chose  1 9/  =  90°  and  (pi  =  225°,  i.  .e.  the  incident 
wavevector  forms  an  angle  of  45°  with  the  normal  to  the  surface.  The  polarization 
both  of  the  incident  and  of  the  scattered  field  was  chosen  to  lie  either  along  the 
meridians  (^-polarization)  or  along  the  parallels  ((^polarization)  of  the  large  sphere 
sketched  above. 

In  the  upper  part  of  the  figure  below  we  report,  for  purposes  of  comparison, 
the  patterns  generated  by  single  spheres  of  radius  p  =  450  nm  and  refractive  index 
n  =  2.0  deposited  on  a  plane  substrate  of  Si  of  refractive  index  n"  =  3.85  +  i0.018; 
the  spheres  are  assumed  to  be  embedded  in  vacuo  (n7  =  1).  In  the  lower  part  of  the 
figure  we  report  the  scattered  intensity  from  a  dispersion  of  the  binary  aggregates 
of  the  spheres  described  above  with  random  distribution  of  their  orientations.  In 


the  left  part  of  the  figure  are  shown,  in  /im  ,  the  quantities  for  the  single 

spheres  and  r2(/w)//o  for  the  aggregates,  where  70  is  the  intensity  of  the  incident 
field  and  the  angular  bracketts  denote  orientational  average.  In  the  right  part  of  the 
figure  the  reported  quantities  are  r2/^//o  for  the  single  spheres  and  r2(/v?tj)//o  for 
the  aggregates  as  an  example  of  cross-polarization. 


At  first  sigth  the  patterns  of  the  co-polarized  intensity  from  single  spheres  and 
from  orientationally  averaged  aggregates  look  rather  similar.  Actually,  the  intensity 
at  the  main  peak  is  larger  for  the  case  of  the  aggregates,  as  it  was  expected  because 
an  aggregate  contains  twice  the  refractive  material  of  a  single  sphere.  The  main 
peak  occurs  in  the  direction  of  reflection  both  for  single  and  for  aggregated  spheres 
(t?s  =  90°,  < ps  =  135°.  Nevertheless  the  width  of  the  peak  from  the  aggregates  is 
smaller  that  for  the  single  spheres.  This  suggests  that  the  effect  of  the  multiple 
scattering  processes  that  occur  among  the  aggregated  spheres  are  not  suppressed 
by  the  averaging  procedure.  In  any  case  a  non-vanishing  field  propagates  along  the 
surface,  although  in  the  ca se  of  the  aggregates,  its  intensity  is  lower  that  for  single 
spheres.  As  for  the  cross-polarized  intensity,  we  notice  that  the  pattern  from  single 
spheres  appears  sharper  than  the  one  from  aggregates.  This  effect  that  we  remarked 
even  in  the  case  of  a  perfectly  reflecting  surface,8  can  be  attributed  to  the  averaging 
procedure. 
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Optical  scattering  methods,  based  on  the  differential  scattered  light  measurements,  are 
extremely  useful  in  particle  sizing  and  characterization  and,  despite  a  large  widespread 
of  proposed  experimental  setups,  there  still  exists  a  possibility  to  modify  and  improve 
them  [1-4].  In  this  work  we  present  an  extension  of  a  previously  reported  differential- 
light  scattering  photometer  suitable  for  measurements  over  a  wide  angular  range  of 
scattered  -  intensity  coming  from  either  aqueous  solutions  of  spherical  particles  or 
particles  deposited  over  a  flat  surface  [5,6].  In  this  scheme,  sketched  in  fig.  1,  the  light, 
coming  from  a  linearly  polarized  10  mW  632.8  nm  He-Ne  laser,  is  focused  on  the 
sample  and  the  scattered  portion  is  collected  by  an  ellipsoidal  mirror  (Melles-Griot  02 
REM  11),  selected  by  a  circularly  shaped  mask  ‘and  by  a  diaphragm  and,  finally, 
brought  directly  to  the  CCD  sensor  of  a  common  8-bit  black  and  white  video  camera. 
The  scattered  intensity  is  acquired  as  an  8-bit  image  through  a  PC  acquisition  card 
using  DMA  transfer  over  the  PCI  bus.  The  image  has  a  circular  ring  shape  and  the 
angular  intensity  distribution  is  represented  by  the  pixel  intensity.  In  order  to  enlarge 
the  measurement  dynamic  range,  which  would  be  intrinsically  restricted  to  8-bit,  we 
used  an  I/O  digital  PC  card  to  change  automatically  and  in  rapid  succession  the  camera 
shutter  opening  times  from  1/60  down  to  1/10000  sec.  In  this  way  a  nearly  15-bit 
dynamic  range  can  be  sampled. 
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Fig.  1  -  Schematic  drawing  of  the  experimental  setup. 

The-present  experimental  setup  offers  the  advantage,  over  the  older  one,  of  a  much 
simpler  configuration,  having  eliminated  two  optical  components  (i.e.  two  focalization 
lenses),  and  allows  a  wider  scattering  angle  interval  to  be  investigated.  On  the  other 

A 

hand  its  application  is  limited  to  static  scattering  experiments.  N 


Fig.  2  -  Measured  scattered  pattern  for  several  particle  diameters  (curves  are  scaled  for  sake  of  clarity). 


In  order  to  test  the  apparatus  effectiveness  we  carried  out  measurements  on  high  purity, 
liquid-chromatography  grade  water  solutions  of  polystyrene  spheres  ranging  from  105 


to  993  nm.  The  results,  referred  to  a  laser  polarization  perpendicular  to  the  scattering 
plane,  are  reported  in  fig.  2.  They  show,  for  all  sizes,  an  excellent  agreement  with  both 
the  previously  measured  ones  and  with  the  theoretical  predictions  which  include 
geometrical  corrections  for  both  reflection  and  refraction  processes  [6,7]. 

This  experimental  setup  is  also  suited  for  measurements  of  light  scattered  from 
particles  on  flat  surfaces.  In  fact  we  deposited  from  a  diluted  solution*  993  nm  latex 
spheres  on  the  specular  surface  of  a  conventional  electronic  grade  silicon  wafer.  We 
checked  the  spheres  distribution  by  optical  microscopy,  which  excluded  the  presence 
of  clusters  and  we  estimated  a  number  of  nearly  40  spheres  contained  in  the  scattering 
area  (a  circle  of '90  |im  diameter).  The  preliminary  Jesuits  obtained  for  both 
perpendicular  and  parallel  laser  polarizations  with  respect  to  the  scattering  plane  are 
reported  in  fig.  3. 


Fig.  3  -  Measured  (points)  and  theoretical  scattered  pattern  for  993  nm  sphere  over  a  silicon  surface  for 
both  perpendicular  and  parallel  polarizations,  with  respect  to  the  scattering  plane,  of  the  normal  incident 
laser  beam.  Scattering  angles  are  measured  from  surface  normal. 
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For  comparison  we  report  also  the  theoretical  patterns  calculated  on  the  basis  of  the 
electromagnetic  field  expansion  into  a  series  of  spherical  multipole  fields  [8].  As  is 
clearly  evident  there  is  a  poor  agreement  in  the  details,  while  the  overall  angular 
•scattered  intensity  distribution  seems  to  agree  better.  We  believe  this  is  mainly  due  to 
both  the  preliminary  nature  of  the  experimental  measurements  and  to  necessity  to 
introduce  some  geometrical  corrections  to  the  pure  theoretical  scattered  pattern. 
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1.  Introduction 


In  the  last  years  several  approaches  have  been-  proposed  for  calculating  the  scattering  pattern 
from  surfaces  seeded  by  particles  of  known  morphology  in  view  of  the  interest  of  this  problem  for 
pure  research  as  well  as  for  technological  applications.*  4  For  the  case  of  a  metallic  surface,  that 
can  be  effectively  modeled  as  a  perfectly  reflecting  plane,  image  theory  proved  to  be  a  suitable 
approach.5  In  practice,  one  has  to  solve,  instead  of  the  real  problem,  the  equivalent  problem  of  the 
scattering  from  the  compound  object  that  includes  the  actual  particle  and  its  image  when  they 
are  illuminated  by  the  superposition  of  the  actual  incident  light  and  of  the  light  that  comes  from 
the  image  source.  The  equivalent  problem  turns  out  to  have  computationally  viable  solutions  for 
single  and  aggregated  hemispheres  with  their  flat  face  on  the  reflecting  surface  and  for  single  and 
aggregated  spheres,6,7  as  well  as  for  cylinders  with  their  axis  parallel  to  the  reflecting  surface.8 

Image  theory  is  not  easily  extended  to  the  case  of  a  dielectric  half-space9  except  when  small  par¬ 
ticles  are  considered.3  When  the  small-particle  simplifications  do  not  apply,  the  scattering  problem 
can  be  alternatively  solved  by  direct  imposition  of  the  appropriate  boundary  conditions  both  across 
the  plane  surface' and  "across  the  surface  of  the  particles.  In  this  respect,  the  method  devised  by 
Bobbert  and  Vlieger10  to  deal  with  spheres  deposited  on  a  dielectric  substrate  is  successful  in 
solving  exactly  the  boundary  value  problem  without  invoking  any  limitation  of  principle  on  the 
size  of  the  particles  to  be  dealt  with.  Nevertheless,  these  authors,  in  order  to  get  the  scattered 
field  irr  the  far  zone,  resort  to  an  approximation  that  has  rather  severe  effects  on  the  reliability  of 
their  results:  their  procedure  predicts,  indeed,  that  no  field  propagates  along  the  surface,  contrary 
to  the  results  of  the  experiments  and  of  ab  initio  simulations.11,12 

The  difficulties  of  the  approach  of  Bobbert  and  Vlieger  have  been  overcome  by  the  present 
authors  through  the  formulation  of  a  general  rule  for  the  reflection  of  a  vector  multipole  field  on 
a  plane  surface.13  This  rule,  when  applied  to  the  description  of  the  scattering  from  a  sphere  on 
or  near  a  dielectric  substrate,14  yields  a  conceptually  simple  procedure  to  impose  the  boundary 
conditions  also  across  the  surface  of  the  sphere  and  to  get  for  the  scattered  field  the  correct  behavior 
at  infinity.  The  results  that  stem  from  this  procedure  are  in  agreement  with  the  experiment  and 
with  the  simulations11,12  and  predict,  in  particular,  the  propagation  of  a  non-vanishing  field  along 
the  surface. 

The  purpose  of  the  present  paper  is  to  extend  the  procedure  of  Ref.  14  to  the  calculation  of  the 
scattering  pattern  from  aggregated  spheres  on  or  near  a  dielectric  surface.  The  need  for  such  an 
extension  is  easily  understood  when  one  thinks  that  aggregation  phenomena  often  occur  among 
the  particles  on  a  substrate.  In  fact,  aggregated  spheres  give  origin  to  anisotropic  particles  whose 
scattering  patterns  are  expected  to  be  rather  different  from  the  patterns  from  single  spheres  and 
to  depend  on  the  orientation  of  the  aggregate  with  respect  to  the  incident  field.  In  this  respect  we 
recall  that  the  pattern  of  the  scattered  intensity  from  a  low-density  monolayer  of  identical  particles 
all  oriented  alike  but  randomly  distributed  on  the  surface  is  merely  proportional  to  the  pattern 
froim-a  single  particle;15  the  case  of  a  monolayer  of  particles  whose  orientations  are  randomly 
distributed  around  an  axis  orthogonal  to  the  surface  requires  instead  to  consider  the  average  over 
the  orientational  distribution.  Since  this  averaging  procedure  implies  further  complications  of 
formalism  with  respect  to  the  case  of  fixed  orientation,  in  this  paper  we  deal  with  the  latter  case 
only  and  defer  the  problem  of  the  orientational  average  to  a  further  paper  that  is  presently  in 
progress. 

In  Section  2  we  sketch  the  theory  oi^  which  the  calculation  of  the  field  scattered  by  an  aggregate 
of  spherical  scatterers  on  the  surface  is  based.  Then,  we  discuss  how  the  full  pattern  of  the  scattered 
intensity  can  be  calculated  for  aggregates  of  fixed  orientation.  Further  details  on  the  mathematics 
of  our  procedure  as  well  as  on  the  related  formalism  and  notation  can  be  found  in  Refs.  13  and  14. 

In  Section  3  we  present  our  specific  results  for  binary  aggregates,  the  scattering  from 
whose  component  spheres  was  already  dealt  with  by  other  authors  with  the  help  of  suitable 
approximations.16,17  We  also  present  our  results  for  the  binary  aggregates  of  spheres  whose  ra¬ 
dius  is  chosen  so  that  the  total  volume  of  the  aggregate  equals  the  volume  of  one  of  the  spheres 
referred  to  above. 
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A  few  concluding  remarks  will  be  drawn  in  Section  4. 


2.  Theory 

A.  Statement  of  the  problem 

Let  us  assume  that  a  plane  surface  separates  a  semi-infinite  homogeneous  isotropic  medium  of 
real  refractive  index  n'  from  another  semi-infinite  isotropic  homogeneous  medium  of  (possibly 
complex)  refractive  index  n":  hereafter  we  willTefer  to  the  region  filled  by  the  former  medium  as 
the  accessible  half-space.  We  also  assume  that  a  plane  wave,  of  circular  frequency  w,  propagates 
within  the  accessible  half-space  and  illuminates  a  cluster  of  spheres,  i.  e.  a  group  of  spherical 
scatterers  whose  mutual  distances  are  so  small  that  dependent  scattering  need  to  be  considered. 
In  order  to  calculate  the  scattering  pattern  from  such  an  object  we  need  to  determine  the  total 
field  at  any  point  P  in  the  accessible  half-space:  we  .will  referjo  the  total  field  as  EInt  when  P  lies 
within  any  of  the  spheres  that  compose  the  cluster,  otherwise  we  will  refer  to  the  field  as  EExt. 
Now,  if  no  particle  were  present,  the  total  field,  EEx\  would  be  the  superposition  of  the  incident 
wave-,  EJ,  and  of  the  field  that  has  been  reflected  by  the  surface,  EH;  these  fields  are  related  to 
-  each  other  by  the  reflection  conditions  on  the  plane  surface.  When  a  cluster  is  present  we  need 
to  add  to  EJ  and  ER  the  field  that  is  scattered  by  the  cluster  itself,  E5,  as  well  as  the  field  that, 
after  scattering  by  the  cluster,  is  reflected  by  the  plane  surface,  E^5.  Ultimately 

EExt  -  E7  +  Er  i-  Es  +  EH5,  ■  -  (1) 

and  Es  is  determined  by  imposing  to  EExt  the  appropriate  boundary  conditions  across  the  surface 
of  each  of  the  spheres  in  the  cluster.  In  this  respect  we  recall  that,  on  the  assumption  that  all  the 
fields  depend  on  time  through  the  factor  exp (-iwi),  the  magnetic  field  B  that  is  also  needed  to 
impose  the  boundary  conditions  is  given  by  the  Maxwell  equation 

iB  =  y  V  X  E, 
k 

with  k  =  oj/c,  that  holds  true  both  for  B/nt  and  BExt. 

The  geometry  that  we  adopt  for  our  study  is  depicted  in  Fig.  1.  The  accessible  half-space 
coincides  with  the  region  z  <  0  of  a  cartesian  frame  of  reference,  of  origin  O,  whose  2  axis  is 
characterized  by  the  unit  vector  z:  the  surface  on  which  the  reflection  occurs  is  thus  the  plane 
z  —  0.  We  assume  that  the  center  of  the  spheres,  of  radius  pa  and  refractive  index  na ,  that 
compose  the  aggregate  lie  at  the  points  Ola  whose  vector  position  will  be  denoted  with  R/a.  It 
is  also  convenient  to  define  the  points  0"  that  are  the.images  of  the  centers  of  the  spheres  with 
respect  to  the  surface  and  to  denote  their  vector  positions  as  R".  The  vector  position  of  a  point 
P  will  be:  denoted  as  r  with  respect  to  O,  as  =  r  —  R'a  with  respect  to  0'a  £hd  as  r"  —  r  —  R" 
with  respect  to  0". 


B.  Multipole  expansion  of  the  field 

The  electromagnetic  plane  wave 
E*  =  £0.,uj„  exp(ik/  -r), 

that  propagates  through  the  halfspace  z  <  0,  is  reflected  into  the  plane  wave 
E*  =  EonURr>  exp(ikR  ■  r), 

where  kj  =  n'kkj  and  k r  =  n'A-k/?  are  the  propagation  vectors  of  the  incident  and  of  the  reflected 
plane  wave  and  u jv  and  u^,  are  the  respective  unit  polarization  vectors  whose  index  77  distinguish 
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whether  they  are  parallel  (77  =  1)  or  perpendicular  (77  =  2)  to  the  plane  of  incidence,  i.  e.  to  the 
plane  that  contains  k/,  kR  and  the  2  axis.  Our  actual  choice  for  the  polarization  vectors  is 

u/i  x  U/2  =  k /,  uhi  x  nR2  =  kR 

with  ur2  =  u/2 .  The  reflection  condition  yields  the  relation 

E0r]  =  Frji^d  ^Eorjj 

where  is  the  angle  between  k /  and  z  and  the  quantities  are  the  Fresnel  coefficients18  for 

the  reflection  of  a  plane  wave  with  polarization  along  u/^ 


*!(*/)  = 


v?  cost?/  —  p 
n2  cos  dj  -f-  P  ’ 


cos  $1  —  (3 
cos  7?/  +  /?  ’ 


where  n  =  n"  / nf  and  *  -  -  ,  ,, _ 

P  =  [(n2  —  1)  +  cos2  7?/] 1/F2. 

Accordingly,’ we  rewrite  the  reflected  plane  wave  as  _  . 

E*  =  Frj{'di)EovuRriexp{ikR  -r). 

For  our  purposes  it  is  convenient  to  expand  both  E7  a'nd  E*  in  terms  of  vector  multipole  fields. 
With  reference  to  the  origin  O  such  expansions  read 

K  ~  E°V  Y1  Jlm  (r  >  n'k)Wlm  >  £/) ,  (2) 

p/m 

E*  =  Jlm(p.»,*)W’«m(fi*».fc*).  (3) 

p/m 

where  we  define  the  vector  multipole  fields18 

J^(r,  K)  =  ji(Kr)X.im  (f),  J^(r,  K)  =  ^Vx  j,(A>)X,m(r)  (4) 


and  the  amplitudes13 

W^(u,K)  =  4nip+1~ 1  (— )m+1  z{Pim  (K)  -u. 

In  the  preceding  equations  the  functions  X/m  are  vector  spherical  harmonics,18 

Z^(K)  =  X;m(K),  Z^(K)  =  X/m  (K)  x  K.  (5) 

are  the  transverse  vector  harmonics  and  the  superscript  p  =  1,  2  is  a  parity  index  that  distinguishes 
the  magnetic  multipoles  (p  =  1)  from  the  electric  ones  [p  —  2).  The  reflection  condition  yields  the 
relation7 

wtfliuR^kn)  =  (— )’,+p+,+m  wfyn  (u/lJ ,  k/) , 

so  that  the  amplitudes  of  E^  never  need  to  be  calculated  explicitly. 

Since  E 7  and  E^  must  satisfy  the  boundary  conditions  also  across  the  surface  of  each  of  the 
scattering  spheres,  we  rewrite  their  expansions,  Eqs.  (2)  and  (3),  in  terms  of  spherical  vector 
multipole  fields  with  origin  at  the  center  of  the  a-th  sphere,  0^,  as  - 
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(6) 


=  EorjUir,  exp[j'k/  •  (R'a  +  r'J] 

=  exp(ikj  •  n'k)W^ (nInM 

plm 

E*  =  Fn(i)j)Eor,URrl  exp[jkR  •  (R'a  +  r'J] 

=  exp(ikfi  •  R !a)Fn (tf/)£'o,  ^  J(m  (ra>  n'k)wi£ (u Rt> ,  kn)-  (7) 

plm 

The  field  that  is  scattered  by  an  aggregate  of  spheres  embedded  within  a  homogeneous  medium 
can  always  be  described  as  the  superposition  of  the  fields  that  are  scattered  by  each  one  of  the 
spheres.  Since  the  latter  fields  must  satisfy  the  radiation  condition  at  infinity  we  write 

Hlm  (r'c  n'k)A{nllm .  (8) 

or  plm 

* 

where  the  multipole  fields  H  are  identical  to  the  J  fields,  Eq.  (4),  except  for  the  substitution  of 
the  spherical  Hankel  functions  of  the  first  kind  ht  for  the  Bessel  functions  ji ;  the  amplitudes  A, 
that  are  as  yet  unknown,  are  determined  by  the  boundary  conditions  at- the  surface  of  each  of  the 

component  spheres.19  *  ~  ~ 

We  still  need  to  consider  the  field  E*5  that,  after  scattering  by  the  cluster;  is  reflected  by  the 

surface.  According  to.  our  previous  results, 1^’1^  at  any  point  of  the  accessible  half-space,  E^  can 
be  written  as 

Ef  -  E0v  £  Him  (C 

or  plm 

i.  e.  as  a  superposition  of  H  multipole  fields  with  origin  at  the  image  points  Oa.  The  amplitudes 
A^Jim  are  related  to  the  amplitudes  of  the  scattered  field  E5  by  the  equation 


~A  (P)  _  V  JP,P')  4(P') 

rjalm  /  jaa:l.l,:m‘nnotl,m' 
p'V 

The  explicit  expression  of  the  amplitudes  is 


(10) 


a 


(p,p') 

a;t,l';m 


—  2Li\n«  htlf,;m^oc;l,,tl,)n 


p"l" 


where  %~l  is  the  inverse  to  the  matrix  ?{(Raa,  n'k)  that  effects  the  transfer  of  the  origin  of  the 

H  multipole  fields  from  O"  to  0’a  and  ROQ  =  R'„  -  R'a'.  In  turn,  the  quantities  are  the 

elements  of  the  matrix  Ta  that  effects  the  reflection  on  the  plane  surface  of  the  H  multipole  fields 
with  origin  at  0'a. 13,14 

The  last  expansion  that  we  need  to  consider  it  that  of  the  field  within  the  spheres  of  the  cluster . 
Since  the  field  within  the  a-th  sphere  must  be  regular  at  0'a  its  expansion  can  be  taken  in  the  form 


plm 


(11) 


C.  Amplitudes  of  the  scattered  field 

In  order  to  calculate  the  amplitudes  A^lm  by  imposing  the  boundary  conditions  across  the  surface 
of  each  of  the  spheres  that  compose  the  cluster,  a  few  steps  are  still  necessary.  When  we  come 
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to  consider  the  a-th  sphere,  the  fields  E(J,  Eq.  (6),  and  ,  Eq.  (7),  that  contribute  to  E^  , 
•  Eq.  (1),  as  well  as'E“,  Eq.  (11),  are  all  given  in  terms  of  vector  multipole  fields  with  origin  at 
0'a.  The  scattered  field  Ejf,  Eq.  (8),  and  the  reflected  scattered  field  E*s,  Eq,  (9),  that  are  given 
by  superpositions  of  vector  multipole  fields  with  different  origins,  can  also  be  expressed  in  terms 
of  multipole  fields  with  origin  at  0'a  through  the  use  of  the  appropriate  addition  theorem.20  The 
result  for  the  scattered  field  E„  is13 


p'l'm'  p?a  J 

where  R^a  =  R'a  -  R^,  while  for  the  reflected  scattered  field  ERS  we  obtain 

plm  p'l'm1  . 

+  E  j 

J 

where  R^  =  R'a  —  R^.  Using  the  expression,  Eq.  (10),  for  the  amplitudes  A  and  taking  into 
account  that  n'k)8mm>  we  get  the  final  result 

J?RS  _  rp  V"^  lOO/V  n'LA  V  A  'r(P*P')  A?  ) 

Erj  — •C'Or;  /  j  **  lm  vro:?  n  )  /  j  ^  \m^7)OcV  m 
plm  p‘l‘ 

+  53  E  5 

P^a  pnlumf 


where  we  define 


We  are  now  ready  to  impose  the  boundary  conditions  across  the  surface  of  each  of  the  scattering 
spheres  by  applying  to  E  and  B  the  procedure  that  we  described  elsewhere  in  full  detail.19  This 
procedure  yields,  for  each  a,  p,  l  and  m,  four  equations  among  which  the  amplitudes  of  the  internal 
field  C  can  be  easily  eliminated.  As  a  result  we  get  for  the  amplitudes  A^lm  the  system  of  linear 
non-homogeneous  equations 


E  E  =  ->v! 


P  p'l'm1 
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and 


R 


(p)  _ 

al  — 


(1  +  na<5pi)u|(nttfc/)a)m(n/fepct)  -  (1  +  n0^2)^/(na^pa)u|(n'fe/?a) 
(1  +  naSpi)u',(nakpQ)wi(n,kfia)  -  (1  +  na6J>-i)ui{nakpa)w'l{nlkpa) 


na  =  ^~  —  1,  ui(x )  =  xjt(x),  wi(x)  =  xhi(x). 
n 

The  quantities  and  coincide  with  the  Mie  coefficients  bi  and  a/,  respectively,21  for  a 
homogeneous  sphere  of  refractive  index  nQ  embedded  into  a  homogeneous  medium  of  refractive 
index  nf.  In  turn 

W$m  =  exp(ikj  •  KWtt(ui„kr)  +  F.itij)  exp(ikR  •  R'a)W^(uRn,kR). 


D.  Scattered  intensity _ _ 


Once  the  amplitudes  A^lm  of  have  been  calculated  by  solving  Eq.  (12),  the  reflected-scattered 
field,  E*5,  is  also  determined  by  Eqs.  (9)  and  (10).  Therefore,  the  superposition  of  E^,  Eq.  (8), 
and  of  E^5,  Eq.  (9),  yields  the  field  that  would  be  observed  by  an  optical  instrument  at  any 
point  in  the  accessible  half-space.  Hereafter  we  disregard  the  reflected  field  because  it  can  be 
observed  in  the  direction  of  reflection  only.  Furthermore,  it  is  convenient  to  use  again  the  addition 
theorem  of  ref.  20  to  refer  both  and  E^5  to  a  common  origin  that  we  choose  to  be  the  point 
O  that  is  defined  in  Fig.  1.  Accordingly,  the  field  that  would  be  observed  in  the  far  zone  of  the 
accessible  half-space  takes  on  the  form 

Ef”  =£0,EH<m(r,riaKlm. 

plm 

with 

At  =  E  E 

a  p'l’m' 

where  the  quantities  J  are  the  elements  of  the  matrix  that  effects  the  transfer  of  origin  of  the  J 
multipole  fields.13  Then,  on  account  of  the  asymptotic  form  of  the  H  fields  for  large  values  of  n'kr 
and  of  their  transversality  in  the  far  zone,  we  are  led  to  write  E^65  as 


j(p') 


(p.p'l 


^(p') 


(13) 


t^06s  exp  (in'kr)^  ^ 

&r)  . —  “  -C'OrjTr?) 


where  we  define  the  normalized  scattering  amplitude18  that,  in  terms  of  the  transverse  vector 
harmonics  Z^(r),  Eq.  (5),  reads14 


H7+'z,c?>( 


(P) 
rjlm  ' 


Therefore,  the  intensity  that  would  be  detected  in  the  direction  k o  with  polarization  along  uor?'  is 

=  ^2  =  ^2  for)  j 

where  Iqv  =  | ^or?  |2  and 
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=«,  -<W  = 


plm 


(14) 


The  efficiency  of  Eq.  (14)  can  be  improved  when,  as  we  assumed  throughout,  the  refractive  index 
v!  is  real.  In  this  case,  indeed,  a  consequence  of  the  addition  theorem  of  ref.  20  is  the  identity13 

exp(-ik-R)WI^(uJk)=  Y,  2M'*).  (15) 

p'/'m' 


where  k  =  n'k k.  Now,  with  the  help  of  Eq.  (15),  the  multipole  transfer  elements  J  that  are 
included  into  the  definition  of  the  amplitudes  A,  Eq.  (13),  can  be  eliminated,  so  that  Eq.  (14)  can 
be  rewritten  as 

U'r,  =  " ixrfk  £ ^  W‘™r (Qorj' ’ ^ 

a  plm  • 

[exp(-ik0  •  Ra)i4^,m  +  exp(-ik0  •  R")^pJ,m] .  (16) 

Equation  (16)  is  the  one  that  we  actually  used  in  our  calculations. 


3.  Results  and  discussion 

The  theory  of  Section  2  has  been  applied  to  the  calculation  of  the  full  pattern  of  the  scattered 
intensity  from  a  binary  aggregate  of  identical  contacting  spheres  embedded  in  vacuo  (n'  =  1)  and 
deposited  on  a  substrate  of  Si  ( nn  =  3.85  +  i  0.018).  We  considered  spheres  of  Si3N4  with  refractive 
index  =  2.0  and  radius  pi  =  450.0  nm,  illuminated  by  radiation  of  wavelength  A  =  632.8  nm. 
The  scattered  intensity  from  single  spheres  of  such  a  morphology  deposited  on  the  same  substrate 
has  already  been  calculated  by  Taubenblatt  and  Tran  through  the  coupled-dipole  method16  and  by 
Johnson17  with  the  help  of  the  Yousif-Videen  approximation.4,22  Our  aim  is  thus  to  use  our  exact 
procedure  to  investigate  the  effect  of  the  aggregation  on  scatterers  whose  properties  are  already 
well  known.  We  also  dealt  with  the  pattern  from  binary  aggregates  of  contacting  spheres  of  SisN4 
with  radius  p2  =  357.0  nm  deposited  on  the  surface  because  the  total  volume  of  such  an  aggregate 
equals  that  of  a  single  sphere  of  radius  p\.  This  choice  allows  us  to  investigate  also  the  effect  of 
the  subdivision  of  the  material  of  a  Sphere.  Hereafter  we  will  refer  to  the  aggregate  of  spheres  of 
radius  pi  as  the  case  of  aggregation  and  to  the  aggregate  of  spheres  of  radius  p2  as  the  case  of 
subdivision. 

The  geometry  that  we  adopted  to  display  our  results  is  the  one  depicted  in  Fig.  1  of  Ref.  6.  Thus, 
the  surface  of  the  substrate  coincides  with  the  zx  plane  and  the  normal  to  the  surface  coincides 
with  the  y  axis.  When  the  single  sphere  is  considered  its  center  lies  on  the  y  axis  while  when  the 
aggregate  is  considered  the  contacting  point  of  the  two  spheres  lies  on  the  y  axis  and  the  axis  of 
the  aggregate  is  taken  to  be  parallel  either  to  the  z  axis  or  to  the  x  axis.  Nevertheless,  the  patterns 
that  we  actually  report  in  this  paper  refer  to  the  former  case  only  because  the  patterns  for  the 
latter  case  do  not  present  specific  features  that  require  a  separate  comment.  The  plane  of  incidence 
coincides  with  the  xy  plane  and  the  direction  of  incidence  forms  an  angle  of  45°  with  the  normal 
to  the  surface;  more  precisely,  the  polar  angles  of  the  direction  of  incidence  are  t ?/  =  90°  and 
Kpi~  225°.  The  incident  light  is  polarized  either  along  the  meridians  (^-polarization)  or  along  the 
parallels  (^polarization)  and  both  the  ^-polarized  and  the  ^-polarized  component  of  the  scattered 
wave  are  considered;  in  practice  the  quantity  that  we  calculated  is  —  r2/^^//or?  and  both  co¬ 
polarized  (lyy  and  1m)  and  cross-polarized  (1##  and  1 vt?)  patterns  are  reported.  The  numerical 
results  from  which  all  the  patterns  were  drawn  achieve  a  convergency  to  four  significant  digits;  this 
required  to  extend  the  multipole  expansions  in  Sections  2  and  3  up  to  l  =■  12.  We  also  recall  that 
all  the  patterns  that  we  are  going  to  discuss  have  a  number  of  features  that  do  not  depend  on  the 
system,  under  study-  but  rather  stem  from  the  boundary  conditions  and  from  the  geometry  of  the 
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problem.  The  occurrence  of  these  features  that  were  discussed  in  Refs.  6  and  14  should  be  borne  in 
mind  for  an  appropriate  interpretation  of  our  results.  Here  we  only  recall  that,  since  we  cpnsider 
the  scattered  field  also  out  of  the  plane  of  incidence,  cross-polarization  effects  are  to  be  expected.6 

In  Fig.  2  we  report  our  results  for  the  pattern  from  a  single  sphere  of  radius  p\  in  contact  with  the 
substrate.  The  section  of  the  X^  and  of  the  patterns  along  the  plane  i9j  =  90°  -can  be  directly 
compared  with  the  results  of  Taubenblatt  and  Tran16  and  of  Johnson.17  The  differences  that  emerge 
from  this  comparison  are  easily  explained  along  the  lines  of  Ref.  14,  where  we  stressed  how  the 
approximations  used  by  other  authors,  possibly  in  conjunction  with  the  approximation  used  by 
Bobbert  and  Vlieger10  to  get  the  far  zone  field,  severely  affect  the  reliability  of  the  calculations. 
The  results  in  Fig.  3  refer  to  a  binary  aggregate  of  spheres  of  radius  pi\  the  patterns  from  the 
aggregate  of  two  spheres  of  radius  p2  are  presented  in  Fig.  4. 

We  first  notice  that  in  agreement  with  our  previous  findings,  the  X^  patterns  show  a  non¬ 
vanishing  field  that  propagates  along  the  surface.  Of  course,  its  intensity  is  smaller  than  the 
intensity  that  we  reported  in  our  previous  paper14  but  this  is  due  to  the  fact  that  the  refractive 
index  of  the  substrate  has  a  non-vanishing  imaginary  part. 

To  confirm  this  interpretation  we  compare  the  pattern  of  from  a  single.sph.ere  in  Fig.  2a  with 
the  pattern  from  the  same  sphere  on  a  perfectly  reflecting  surface.  According  to  our  calculations 
the  latter  pattern  (not  reported  here)  has  a  rounded  shape  in  contrast  to  the  rather  peaked  shape  of 
the  pattern  in  Fig. '2a:  We  also  remark  that  also  the  patterns  that  we  reported  in  Fig.  5  of  Ref.  14 
for  spheres  of  different  size^and  refractive, index  and.  for  several  choices  of  the  (real)  refractive 
index  of  the  substrate  present  a  rounded  shape-.  We  are  thus  led  to  conclude  that  a  non-vanishing 
imaginary  part  of  n"  produces  a  strong  damping  of  the  scattering  patterns.  As  a  result  we  are  left 
with  a  noticeable  peak  around  the  direction  of  reflection  {d$  =  90°,  <ps  —  135°).  This  effect  is 
not  restricted  to  the  case  of  a  single  sphere  because  the  same  peaked  shape  of  the  X and  of  X$$ 
patterns  can  be  observed  in  the  case  of  the  aggregated  spheres;  in  other  words  this  feature  does 
not  depend  on  the  shape  of  the  scatterers  and  should  thus  be  due  to  the  non-vanishing  imaginary 
part  of  n" .  The  height  of  the  main  peaks,  for  ^-polarized  incident  field,  depends  more  or  less  on 
the  quantity  of  refractive  material:  for  instance,  the  height  of  the  X peak  is  almost  equal  for  the 
single  sphere  and  for  the  case  of  subdivision  whereas  the  height  for  the  case  of  aggregation  is  about 
four  times  larger.  Furthermore,  the  X ^  pattern  for  the  case  of  aggregation  show'a  noticeable 
increase,  with  respect  to  the  case  of  subdivision,  in  the  direction  ds  —  90°,  (ps  =  180°. 

Anyway,  the  largest  value  of  the  scattered  intensity  does  not  occur  exactly  in  the  direction  of 
reflection  (tf5  =  90°,  (ps  =  135°).  This  may  be  surprising  because,  according  to  our  calculations 
for  the  same  particles  in  free  space,  the  largest  scattered  intensity  occurs  in  the  forward  direction 
that,  for  particles  on  a  surface,  coincides  with  the  direction  of  reflection.7  Actually,  for  the  single 
sphere  the  maximum  both  of  Xvv  and  X **  occurs  at  <ps  =  155°,  so  that  the  coupling  with  the 
surface  is  so  large  as  to  produce  a  significant  shift  of  the  main  peak  with  respect  to  the  direction  of 
reflection.  For  the  aggregate  of  spheres  of  radius  p\  the  maximum  of  Xw  occurs  at  (ps  =  130°  but 
it  occurs  at  ips  -  150°  for  I**;  finally  for  the  aggregate  of  spheres  of  radius  p2  the  maximum  of 
X^y  occurs  again  at  (ps  —  135°  but  the  maximum  of  X$$  occurs  skips  =  150°.  Even  though  we  do 
not  report  the  specific  patterns  we  can  state  that  a  similar  behavior  occurs  even  when  the  axis  of 
the  aggregates  is  parallel  to  the  x  axis.  This  seemingly  erratic  location  of  the  maximum  arises  from 
the  coupling  of  the  scattering  particles  with  the  substrate  and,  for  the  case  of  aggregated  spheres, 
also  on  the  competing  coupling  of  the  component  spheres  among  themselves.  The  strength  of 
both  these  couplings  is  likely  to  differ  for  different  polarization  whereas  the  coupling  among  the 
component  spheres  is  expected  to  depend  on  their  size.  Hence,  the  dependence  of  the  location 
of  the  maximum  on  the  polarization  as  well  as  on  the  radius  of  the  aggregated  spheres  can  be 
explained  in  terms  of  the  ratio  between  these  couplings. 

Finally  we  notice  that  the  XV(p  and  1#$  patterns  from  single  spheres  as  well  as  from  aggregated 
spheres  both  of  radius  pi  and  p2  are  of  rather  similar  shape,  even  though  the  specific  values  of  the 
scattered  intensity  may  be  noticeably  different  as  we  noted  above  with  regard  to  the  maximum. 
This  means  that  an  analysis  of  the  scattered  intensity  in  the  plane  of  incidence  only  cannot  yield 
enough  information  to  discriminate  the  shape  and  size  of  the  particles.  In  turn,  even  a  cursory 
examination  of  the  cross-polarized  patterns  X ^  and  X^$  show  that  they  depend  both  on  the  shape 
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and  on  the  size  of  the  scattering  particles.  We  are  thus  led  to  conclude  that  any  attempt  to  gain 
information  on  the  morphology  of  particles  deposited  on  a  dielectric  substrate  requires  considering 
both  the  co-polarized  and  the  cross-polarized  aspects  of  their  scattering  pattern. 


4.  Conclusions 

The  patterns  that  we  discussed  in  Section  3  are  only  a  sample  of  the  changes  undergone  by  the 
scattered  intensity  when  spheres  deposited  on  a  dielectric  substrate  aggregate.  However,  the  results 
that  we  reported  prove  that  the  present  extension  of  the  theory  of  Ref.  14  to  aggregated  spheres 
preserves  the  distinctive  feature  of  the  original  approach  that  a  non- vanishing  field  can  propagate 
along  the  interface.  We  stress  again  that  this  feature  is  due  to  the  fact  that  our  procedure  is  based 
on  no  approximation,  neither  to  get  the  far  zone  expression  of  the  scattered  field.  Strictly  speaking, 
our  approach  requires  the  truncation  of  the  multipole  expansions  on  which  the  formalism  is  based. 
Nevertheless,  this  unavoidable  truncation  is  of  numerical  and  not  of  physical  nature.  Therefore, 
once  the  convergency  has  been  checked,  our  results  remain  physically  reliable. 

A  noticeable  advantage  of  any  approach  based  on  the  multipole  expansion  of  the  field  is  the  pos¬ 
sibility  of  using  the  transformation  proper ties,  of  the  ^vector  multipole  fields  under  rotation.23  Once_ 
the  scattering  amplitude  of  an  anisotropic-scatterer  has  been  calculated  for  a  given  orientation, 
producing  the  pattern  for  any  other  orientation  is  a  fast  and  low-cost  operation  thaL  produces  a 
large  amount  of  data.1 2 3 4 5 6 7 8 9 10 11 12  Nevertheless,  in  the  present  paper  we  were  interested  in  the  effects  that  are 
due  to  the  anisotropy  of  the  aggregates  so  that  we  resolved  to  present  the  patterns  for  a  specific 
choice  of  the  orientation  and  of  the  direction  of  incidence,  and  to  defer  to  a  forthcoming  paper  the 
investigation  of  the  patterns  from  randomly  oriented  aggregates. 
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Fig.  1.  Sketch  of  the  geometry  that  we  adopted  in  our  theory.  Only  the  a-th  sphere  of  the  cluster  is 

shown  for  the  sake  of  clarity. 


Fig.  2.  Pattern  of  the  scattered  intensity  from  the  single  sphere  of  radius  pi.  The  quantity  that  is 
actually  reported  (in  pm2)  is  T<p<p  in  (a),  X-q^  in  (b),  T ^  in  (c)  and  Z00  in  (d). 


Fig.  3.  Pattern  of  the  scattered  intensity  from  the  binary  aggregate  of  the  spheres  of  radius  pi.  The 
quantity  that  is  actually  reported  (in  pm2)  is  (a)>  in  (b)»  (c)  an<^  ^-d  in  (d). 


Fig.  4.  Pattern  of  the  scattered  intensity  from  the  binary  aggregate  of  the  spheres  of  radius  p2.  The 
quantity  that  is  actually  reported  (in  pm2)  is  X<p<p  in  (a),  Xo^,  in  (b),  X in  (c)  and  X-o-o  in  (d). 
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